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Abstract 

The purpose of this paper is to establish LP error estimates, a Bernstein in- 
equality, and inverse theorems for approximation by a space comprising spheri- 
cal basis functions located at scattered sites on the unit n-sphere. In particular, 
the Bernstein inequality estimates LP Bessel-potential Sobolev norms of func- 
tions in this space in terms of the minimal separation and the LP norm of the 
function itself. An important step in its proof involves measuring the LP sta- 
bility of functions in the approximating space in terms of the £ p norm of the 
coefficients involved. As an application of the Bernstein inequality, we derive 
inverse theorems for SBF approximation in the L p norm. Finally, we give a new 
characterization of Besov spaces on the n-sphere in terms of spaces of SBFs. 



1 Introduction 

Various applications in meteorology, cosmology, and geophysics require a modeling of 
functions based on scattered data collected on (or near) a sphere; i.e., when one docs 
not have any control on where the data sites are located pj [SJ On the unit 
sphere in R n+1 , n > 1, a popular method is to construct the required approximation 
from spaces of spherical basis functions (SBFs), which are kernels located at points 
in a discrete set X = {£j}|Li € §>", the set of centers or nodes. 

A function (j> : [— 1, 1] — * R is an SBF on §" if, in its expansion in ultraspherical 
polynomials P( , A„ = the Fouricr-Legendre coefficients {4>(i)} of <f) arc ai l 

positive; see section [3] for details. These 4> are to be used as kernels of the form 
(j)(x ■ y), x, y G S™, x- y being the usual "dot" product. The approximation space here 
is the span 

Q<f,,x ■= span{0(a- • £)} 5e x- 
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Following usage common in the neural network community, we will say that a function 
.9 £ G<t>,x is an SBF network associated with (j>. The SBF cf> is sometimes called an 
activation function or a neuron, but we will not use these terms here. 

Such <j> may have singular behavior. This is the case for certain thin-plate splines; 
(1 — x-y)~ x / 2 is an SBF in n > 2, for instance. However, when they are continuous, 
they are positive definite in Schoenbcrg's sense [30 . In that case the interpolation 
matrix • is positive definite, and it is possible to use SBFs to interpolate 
data given at the points in X . 

The focus of this paper is approximation. To handle noisy data, both least squares 
and quasi-interpolants have been used for many years. More recently, the issue in 
many meshless numerical methods for solving PDEs is how well a network approx- 
imates a solution to the PDE. Singular SBFs should prove useful in probing for a 
corresponding singularity in solutions. 

To be effective, though, such methods require knowing the degree of approxima- 
tion in various spaces, especially the L p , 1 < p < oo. The L 2 case for SBFs with 
4>(£) ~ (£+ 1) -/3 , (3 > n/2 has recently been investigated in [53], with nearly optimal 
rates being attained by interpolatory networks. The known estimates on the degree 
of approximation in the case of L p , p ^ 2 provided by interpolatory networks are 
not asymptotically optimal. This has lead to the development of other approxima- 
tion tools Q21 EH [H] , involving SBFs or spherical harmonics, in L p , 1 < p < oo. 
A central step in obtaining approximation rates in 1? was establishing a Bernstein 
estimate, which was the used to get an inverse approximation theorem. 

The paper has three main goals. The first is to derive an LP Bernstein inequality, 
for 1 < p < oo; namely, ||g||#p < Cq~ 1 \\g\\ p , < 7 < c^. Here H p is a Bessel- 
potential Sobolev space [32l[34j; it measures derivatives of g (cf. section |2~3|) . The 
quantity q is a half the minimal separation of points in X; q^ 1 plays the role of a 
Nyquist frequency. 

The second is to obtain is to obtain L p error estimates, 1 < p < 00, for approx- 
imating a function by networks in G<j,,x- We combine these direct (Favard- Jackson) 
estimates with the Bernstein inequalities to provide new characterizations of Besov 
spaces on S™, characterizations that use rates of approximation from the G^.x- The 
Bernstein estimates are then used to establish inverse theorems and obtain nearly 
optimal rates of approximation. 

The third is to show that the results gotten here will apply for nearly all of the 
SBFs of interest. In particular, they apply to various RBFs restricted to the sphere 
- the thin-plate splines and Wendland functions, whose Fourier-Legendre coefficients 
have algebraic decay, and also Gaussians and multiquadrics, whose coefficients decay 
faster than algebraically. SBFs in the latter class are well known to be difficult to 
treat. 

The paper is organized this way. Section [5] reviews various geometric quantities, 
such as the set of centers, mesh norm, and so on. It also discusses spherical harmonics 
and the Bessel-potential Sobolev spaces. Section [3J discusses SBFs, their Fourier- 
Legendre expansions, and deals in detail with the SBFs mentioned earlier, along 
with ones corresponding to certain Green's functions that play a significant role in 
the paper. It is here that we will show that nearly all of the SBFs of interest have 
the properties necessary for our results will hold. We also mention that we obtain 
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precise asymptotic expressions for the Fourier-Legendre coefficients in the case of the 
Wendland functions. 

The strategy for establishing the Bernstein inequality, which will be detailed be- 
low, consists of two key components: LP approximation results for functions in Q^.x 
by means of spherical polynomials, and LP stability estimates; these are developed in 
scctions[2]and[5J respectively. The approximation results are based on Marcinkicwicz- 
Zygmund inequalities developed in [17l [16j [21] , as we U as frame results from J5JJ. The 
stability results, which are of interest in their own right, are for all LP - not just for 
interpolation with continuous SBFs. To obtain them, we introduce a stability ratio, 
which provides some measure of the extent to which a finite set in LP is linearly 
independent. 

In section [6j the results of the previous two sections are combined to yield LP 
Bernstein inequalities f section 16. 1| . direct theorems for approximation by networks 
in G<p,x (section ROj) , characterizations of Besov spaces on S™ (section RTB"! , and inverse 
theorems for LP functions approximated at given rates by SBF networks (section l6.4l) . 

Strategy Let g be an SBF network in Gc/>,x C H^(S n ), so that it has the form 

#( x ) = H M( x "0- 

One of our main goals is to obtain an LP Bernstein inequality for such networks; that 
is, a bound of the form ||<7||#p < C<7~ 7 || <?|| p , where the norms are those appropriate 
for S™ and 7 > is bounded above by a constant depending on <fr and p. 

Our strategy involves approximating g by degree L spherical polynomials on 
where L ~ q . Now, for fixed L and any S, there is a Bernstein inequality, ||S||#p < 
CL 1 1| S\\ p , which is found in Theorem 14. 101 Using it and manipulations involving the 
triangle inequality, one has that 

< + 11.9 - s\\ H » < CL-y\\s\\ p + \\g- S\\ H *, 

which holds for given L and any S. 

Obtaining an appropriate polynomial S is crucial to the argument. To do that, 
we will use the frame operators introduced in [21| and discussed in more detail in 
section 14.31 below. In particular, we need reconstruction operators Bj, with J ~ 
log 2 L. These rotationally invariant operators have other very useful approximation 
properties, which arc given in Proposition 14.91 They take LP spaces and the space 
of continuous function boundedly into spherical polynomials having degree 0(2 J ). 
Consequently, with S = Bjg, we have \\S\\ P < C||c/|| p , and also 

IMk ? < C2^\\g\\ p + ||.9 - Bjg\\ HPy = CT J \\g\\ p + jj- • l|g ~ ^ g|l ^ P . || ff || p 

ll.9llp \ a \p 

where \a\ p = ( z^^ex \ a i\ P j i s the p-norm of a = {a^}^x- 

The functions {</>((•) ■ are linearly independent and form a basis for Q, and 

so the pairing a <-> g is bijective. Since Q has finite dimension |^|, the ratio 

rg iP := max — -f- (1.1) 

GBg^O g p 
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is finite; it will be called the p-norm stability ratio of the network Q = Q$,x- This ratio 
is similar to a condition number in interpolation, but for L p . With it, the inequality 
directly above becomes 

IMk < (C2i j + C'r g , p ( " (J ~|^ )g "^ )) \\9\\v (1-2) 

To obtain the desired Bernstein inequality, we require two bounds: the first on 
||(/ — B ,7)5!! H p/\a\ p and the second on fg^ p . The first bound relies only on approxi- 
mation results; these we cover in section |4l The second is a bound on the stability 
ratio. This bound requires a more detailed analysis involving both the geometry of 
X and properties of </>. It is carried out in section [5] 

An interesting point is the that the two bounds make different demands on the 
properties required for 0. This makes the analysis of both bounds subtle. Fortunately, 
the common demands are satisfied by large classes of SBBs, including restrictions to 
S™ of the most common RBFs - the thin-plate splines, Wendland functions, Gaus- 
sians, Hardy multiquadrics, and others. 

2 Background 

2.1 Background and notation for S n 

Centers and decompositions of S n . Let X be a finite set of distinct points in 
we will call these the centers. For X, we define these quantities: mesh norm, 
hx = sup^ggn inf j 6j jf g?(£, y), where d(-,-) is the geodesic distance between points 
on the sphere; the separation radius, qx = \ min^^/ £'); and the mesh ratio, 
Px ■= hx/qx > 1- 

For p > 1, define T p = JF P (S") be the family of all sets of centers X with px < p ■ 
We say that X is p-uniform if X G T p . For every p > 2, F p (E> n ) is not only non empty, 
but it contains nested sequences of sets of centers for which hx becomes arbitrarily 
small; precisely, the result is this: 

Proposition 2.1 ( |23L Proposition 2.1]) Let p > 2 and let T p be the correspond- 
ing p-uniform family. Then, there exists a sequence of sets Xk € T p , k = 0, 1, . . 
such that the sequence is nested, Xk C Xk+i, and such that at each step the mesh 
norms satisfy \hx k < hx k+1 < \hx k - 

We will need to consider a decomposition of §™ into a finite number of non- 
overlapping, connected regions R^, each containing an interior point £ that will serve 
for function evaluations as well as labeling. For example, if X is the Voronoi tes- 
sellation for a set of centers X, then we may take to be the region associated 
with £ G X. In any case, we will let X be the set of the £'s used for labels and 
X = {R 6 C S n I £ G X}. In addition, let = max CeA -{diam(i? c )}. 
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2.2 Spherical harmonics 

Let n > 2. Let d/x be the standard measure on the n-sphere, and let the spaces L p (§"), 
1 < P < oo, have their usual meanings. In addition, let Ag™ denote the Laplace- 

Beltrami operator on S™. The eigenvalues of Ag™ are —£(£ + n — 1), £ = 0, 1, For 

n > 2 and fixed, the dimension of the eigenspace is 

]n _ £ + \ n fl + n-2\ T- 1 _n-l 

d< "~\, / J ~ A„(n-2)!' ' _ 2 (2 ' 1} 

For n = 1, the case of the circle, dp = 1 and d] = 2, £ > 1. 

A spherical harmonic Y^ m is an eigenfunction of Ag™ corresponding to the eigen- 
value -£(£ + n - 1) Pi, where m = 1 . . . <%. The set {Y Ljn : I = 0, 1, . . . , m = 
1 . . . d™} is orthonormal in L 2 (S n ). Denote by Tig the span of the spherical harmonics 
with fixed order £, and let II l = ®£_ Tit be the span of all spherical harmonics of 
order at most L. The orthogonal projection Pi onto Tig is given by 

Ptf= X)</>*WVi,™- (2.2) 

m— 1 

We regard the sphere §™ as being the unit sphere in R™ +1 , and we let the quantity 
£•77 denote the usual "dot" product for R" +1 . Using the addition formula for spherical 
harmonics, when n > 2, one can write the kernel for this projection as 

m ■n) = Y J YeMOY^n) = £ -±±ip^\t ■ V ), x n ■.= (2.3) 

m— 1 

where Pi (•) is the ultraspherical polynomial of order A„ and degree £. Also, we 

have that H-F^ ||oo < F^ (1) = fe^ 2 ■ Wc will briefly discuss these polynomials in 
section [31 in connection spherical basis functions. For n = 1, Ai = 0. In that case, 
the kernel for Pi has the form 

m-v) = \ f . />1 (2-4) 

where Tf(-) the degree-^ Chebyshev polynomial of the first kind, which is a limiting 
case of the ultraspherical polynomials [331 Section 4.7]. 

We will also need to consider operators of the form J2 < eLo c ^ e - The kernels for 
the projections Pi then provide us with kernels X^=o c t^(^ ' v)> which may be dis- 
tributional. 

2.3 Bessel-potential Sobolev spaces 

The spherical harmonic Yg^ m is an eigenfunction corresponding to the eigenvalue 
— £{£ + n — 1) = \\ — {£ + A„) 2 for Laplace-Bcltrami operator Ag™ on S". It follows 



5 



that I + A„ is an eigenvalue corresponding to the cigcnfunctions Yi iTn , m = 1 . . . d™, 
of the pseudo-differential operator 

oo 

U := VK- As- =X> + A «) P *- ( 2 - 5 ) 

Let 7 be real, 1 < p < oo and n > 2. If / is a distribution on §™, define the 
Bessel-potential Sobolev spaces HP(E> n ) [32j[34] to be all / such that 

oo 

ll/llffi? := |X> + A «W r <0 °' ( 2 - 6 ) 

T II ' — ' Lp 

e=o 

where Pg is from (|2.2|) . The notation we use here is that of Triebel [33]. Strichartz 
32] defined these spaces on an a complete Ricmannian manifold, using the equivalent 
operator (1 — Ag™) 7 / 2 to do so. One more thing: 

Remark 2.2 The space iJ 2 (§ Tl ) is the domain o/L 7 [32j Theorem 4.4], which implies 
that .ff 2 (S™) is norm equivalent to the usual sobolev space W7(S n ). 



3 Spherical basis functions 

For any real A > 0, not just A n = ^-j^, the ultraspherical polynomials satisfy the 
orthogonality relation, 

For the circle, we have Ai = 0. With I > 1, as A — > 0, the ratio P^\-)/X converges 
to (2/£)Te(-), the degree-^ Chebyshev polynomial of the first kind [33J Section 4.7]. 

Consider a function <f> in LP or C. We will assume that cf> has the following 
expansion in the orthogonal set of ultraspherical polynomials: 

f 2^(O) + iE^W0, »=1, 
^ I E^Wt^^W), n>2. 

where uo n := prnTn i s t ne volume of S™. 

Functions of this form are called zonal. We will assume that the series converges 
in at least a distributional sense. The coefficients in the expansion are obtained via 
the orthogonality relations in (|3.ip . These are given below. 

i + A n ^ (i + A n )r3(A„)r(i + i) p( A.). V1 



)1-A, 



7rr(£ + 2A„) 



(x)P^" ; (x)(l-a; z ) A "- = cix. 
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Using Rodrigues' formula [33j Eqn. (4.7.12)] for Pg'(x) in the equation above and 
employing the duplication formula and other standard properties of the Gamma func- 
tion, one can obtain this expression: 



dx, 



which holds for all £, even when n = 1 - i.e., Ai = 0. 

Schocnbcrg [3D] defined <fi to be positive definite if for every set of centers X the 
matrix [</>(£, j ■ £,k)] is positive scmidcfmite. He showed that cf> is positive definite if and 
only if the coefficients satisfy (f>(£) > for all £ and J27Lo fii^de < oo. If in addition 
(/>(£) > 0, then [</>(£•/ ■ is a positive definite matrix and one can use shifts of <j) 
to interpolate any function / £ C(S n ) on X. We will say that is a spherical basis 
function (SBF) in this case. 

One usually makes the assumption that the sum X^fco 'fiitydi < °°' f° r then <f> is 
continuous and (f)(1) = ||</>||loo. This is essential if we are doing standard interpolation 
of a function from its values on X. However, we are more interested in approximation 
than interpolation, and so we will not make this assumption here. Indeed, we will 
say that any distribution cf> for which cf>(£) > for all £ is a spherical basis function. 
In general, we will be interested in SBFs in V . 

Zonal functions that satisfy 4>(£) > for £ > L > are said to be conditionally 
positive definite SBFs. In the RBF theory on Euclidean space, the difference between 
strictly positive definite RBFs and conditionally strictly positive definite RBFs is 
significant. On §™, this difference is less important: a conditionally positive definite 
SBF differs from an SBF by a polynomial of degree L — 1. This does play a role in 
interpolation, but is much less significant in approximation problems. That being the 
case, unless there is a genuine need to distinguish between the two, we will refer to 
both as simply SBFs. 

Below we will list Fourier-Lcgcndre expansion coefficients for some of the more 
significant SBFs. Apart from certain Green's functions that we will do first, these 
SBFs are restrictions of Euclidean RBFs in W l+1 to the § n , which are themselves 
SBFs [Ml Corollary 4.3]. These include Gaussians, multiquadrics, thin-plate splines, 
and Wendland functions. Such SBFs are RBFs expressed in terms of the Euclidean 
distance between £ and rj or its square, ||£ — ?7|| 2 = 2 — 2£ • n and, with t = £ ■ ?y, these 
give rise to functions of 1 — t. 

Green's functions Let f3 > 0. The Green's function solution to L^G^ = S is a 
kernel with an expansion in spherical harmonics having coefficients Gp(£,m) = (£-{- 
. Properties of Green's functions are discussed in more detail in Proposition ^. 121 
We simply remark that the kernel Gp is an SBF that is in L 1 (S n ) for all f3 > 0. For 
us, Gp will play a significant role. The SBFs we consider will generally be of two 
types: (j> = Gp + Gp * tp, where ip is an L 1 zonal function, or <\> will be in C°°. The 
first type includes the thin-plate splines and Wendland functions, and the second, the 
Gaussians and multiquadrics. 
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Thin-plate splines The thin-plate splines are defined in [351 Section 8.3]; their 
Fourier-Legendre coefficients are found in [22l §4.2]. These are given below. 



>(*) 



(-ljrw+i^-*)-, 
i-iy+Hi -i)-iog(i-t), 

x (f\-r, T{l - s) 

where the factor C S! „ is given by 

( sill(7Ts) 



s G N. 



(3.3) 



c s ,„ : =2 s +"^r( s + i)r( s + -) 



i. 



s > -§, s g* N 
s G N. 



Let v = £ + X n . For large i/, the Fourier-Legendre coefficients 4> s (£) for the thin- 
plate splines have the asymptotic form 



p-i 



4> s {£) = C StU V 



1 + Gj(n, s)v J +R p (n,s, 

3 = 1 



(3.4) 



where Rp(n, s, v) = 0(z/ p ) and Gj(n,s) are defined in p. 119]. 

Two remarks. First, we have made use of Go(n,s) = 1 in the expansion from 
[25l p. 119]. Second, when s is an integer or half-integer, 4> s (l) is a rational function 
of £, and, hence, of v. In that case, it follows that the series for </) s (£) is actually a 
convergent power series in v . For other s, the expansion is only asymptotic. 

From the structure of the expansions above and the properties of Green's functions 
listed in Proposition 14.121 we see that any finite linear combination of thin-plate 
splines 

m 

4> = y^.Aj(j) Si , -- < si < s 2 < ■ ■ ■ < s m , (3.5) 



has the form 



4> = A!(G 2s+ „ + G 2s+n * V), ^ e L 1 . 



(3.6) 



Wendland functions All of the SBFs we have discussed so far are related to RBFs 
stemming from completely monotonic functions. These RBFs have the property that 
they are strictly positive definite or conditionally positive definite in K™ for all n. 
The corresponding SBFs are also positive definite in again for all n. These RBFs 
are not compactly supported, however. This can be remedied, but there is a price: 
we must give up positive definiteness beyond a certain dimension. 

Wendland (cf. [35l Section 9.4]) constructed families of RBFs that are compactly 
supported on < r < R, strictly positive definite in Euclidean spaces of dimension d 
or less, have smoothness C 2k , and, within their supports, are polynomials of degree 
L^J + 3fc + 1. The quantities d, k, and R are parameters and may be adjusted as 
needed. 
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Restricting the Wendland functions to S n just requires setting r = \/2(l — t) and 
R = \/2(l — to), where — 1 < to < t < 1. We will denote these functions by <frd,k(t). 
The support of </> d , fc on S" is then < 9 < cos" 1 ^) < 7r. Fr om [3 5l Theorems 9.12 
& 9.13], if t > to ; then these functions are polynomials in y/l — t that may be put 
into the form, 

Pd,k(t) = ei(l - t) + (1 - t) fc+ 2e 2 (l - i), 

where ej and e 2 are polynomials having with degei = [^(LfJ +3fc + l)J and dege 2 = 
|_i(|_|J + fc)J ■ Outside of this interval, the <j>d,k are identically 0. Using a power series 
argument, we have that, near t }t to, <j>d,k{t) = A(t - t )^^ +2k+1 (l + 0(t — to)), from 
which it follows that <fid.k{t) is piecewise C^^ +2k+1 near to. In addition, it follows 
that tpd,k{t) ■= <t>d,k{t) —Pd,k{t) i s piecewise C^l+ 2k + 1 on the whole interval [—1, 1]. 
Putting all of this together, we conclude that 

0d,fe(t) = ei(l - t) + (1 - t) fc+ 'e 2 (l - t) + Vd, fe (t). (3.7) 

Our aim is to use this decomposition to obtain large £ asymptotics for the Fourier- 
Legendre coefficients 4>d,k{£) in This we now do. 



Proposition 3.1 Let m = |_f J + 2fc + 1. If £ > degei, then 

4>d,k(i) = {£ + A„)-( 2fc+1+n ) (a + j^- + 0(£ + A„)- 2 ) + 
Moreover, if we choose > n, then the c/>d,k have the structure 



^ d .k{£) 

{£ + Xn) m ' 



<t>d,k = polynomial + A [G 2 k+n + G 2fc+n * tfjj , xj) e L 

Proof: The polynomial term ei(l — t) doesn't contribute to coefficients with £ > 
degei. The term (1 — t) fe+ 2e 2 (l — t) is a linear combination of thin-plate splines, 
starting with s = k + |. Thus it contributes the first term on the right above. By 
Remark [2~2l the function ipd,k is in H^, so it can be written as tpd,k = L„; m L™^,fc- 
The second term on the right follows directly from this fact. Finally, the form of the 
4>d,k(£)'s leads to the second statement. □ 

Before leaving the topic, we point out that, when [_|J > n, we have determined the 
precise asymptotics of the Fourier-Legendre coefficients for the Wendland functions. 
Heretofore only upper and lower bounds were known. 

Gaussians The Fourier-Legendre coefficients for the Gaussians, which arc given 
below, may be found in [Ml Ex. 37, p. 383], QH Example 5.2], and [H §4.3]. 



7ff (t) = e -Mi-t), a >0, 
y„(£) = 2n(^) K e-"h n+e (a), 



(3.8) 



9 



where I\ n +i is an order A„ + £ modified Bessel function of the first kind. For all £ > 0, 
the coefHcient 7<r(£) satisfies this bound: [22j Proposition 4.3]: 



r(i + 2±i) 



r(i + s±i) 



(3.9) 



Multiquadrics The Hardy multiquadrics are treated in [22l §5]. The results are: 



mq a (t) = -^2 + 2 (i-t), tf > 0. 

7r A "r(^- 1/2) 



V + 2)^i/2r(£ + A„ + i) 

7- 1/2 £+1/2 



2^1 



■J + A„ + 1; 



(5 2 + 2)^ 



(3.10) 



Here, 2-F1 is the usual hypergeometric function. Expressions for Fourier-Lcgcndrc 
coefficients for generalized multiquadrics may be found in [221 §5]. Again, this time 
for £ sufficiently large, the coefficient inq s (£) satisfies the following bound [551 Propo- 
sition 5.1]: 



1 



S 2 + 2 



< rhq 5 (£) < C 2 £ 



-l-n 



S 2 + 2 



(3.11) 



Ultraspherical generating functions For n > 2, the ultraspherical polynomials 
Pj A "^ are frequently defined in terms of the generating function [33J Equation (4.7.23)] 
below: 

u\ n , w (t) = (1 - 2tw + w 2 )- Xn , 1 > w > 0, n>2 ) 



U\ n ,w(£) = W £ 



(3.12) 



When n = 1, Ai = 0, the expansion is in terms the Ti(t)'s, the Chebyshev polynomials 
of the first kind. In this case, the gerating function is simply the Poisson kernel. 



P w {t) 



l-2tw+w 

1, * = o, 

2w e , £ > 1 



1 > w > 0, 



(3.13) 



4 Approximation 

The approximation part of the analysis makes use of kernels and frames, which are 
related to them. These were studied in [IJ [T21 HU [H3 [H] and further developed 
in [26j : we review them here, along with a number of other results important to 
attaining the goals of this paper. First, we will develop various types Marcinkiewicz- 
Zygmund inequalities for the sphere. Although some of these were previously derived 
[T71 [TBI |2"T] , those pertinent to both the approximation and stability analysis are new. 
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Second, using frames we establish a Bernstein inequality for spherical polynomials. 
Again, using frames we establish various distance estimates for </> G Hp and we discuss 
Green's function solutions to L^G^ = 5. As we have mentioned earlier, these form a 
very important class of SBFs. Finally, at the end of the section we will complete the 
approximation part of the analysis. 



4.1 Kernels 

Let re(i) G C^R), with k > max{2,n — 1}, be even, not identically 0, and satisfy 

\K {r) (t)\ < C K (1 + \t\) r - a foralUeK, r = 0,...,k, (4.1) 

where a > n + k and C K > are fixed constants. We remark that all compactly 
supported, C k functions that are even satisfy (|4.ip . Functions in the Schwartz-class 
S(M) that are even satisfy (|4.ip for arbitrarily large k and a. Given such a k, define 
the family of operators 

oo 

K £i „ := K (eL n ) = ]T k(e{£ + X n ))P £ , < e < 1, 
along with the associated family of kernels 

j 2¥ K (0) + ^Efcl^(^)cOS^, 71=1, 

KeA iJ" \ J2Zo<e(i + X n ))^P^\cose), n>2, 

COS n n 

where cos 9 = £ ■ r\ and < e < 1 . 

It is worthwhile noting that «(i) = e _t " satisfies (|4.ip and that the corresponding 
kernel is essentially the heat kernel for 

We will need several results concerning these kernels and operators. First of all, 
we require the estimates on the L p norms for the kernels. Material closely connected 
to the theorem below appeared in [TSJ Proposition 4.1]. 

Theorem 4.1 ([HI Theorem 3.5 &; Corollary 3.6]) Let k satisfy (|4~T|) . with k > 
max{2,?7 — 1}. I/O < 9 < tt, then there is a constant (3 n ,k,K > such the kernel K e ^ n 
satisfies the bound 

Pn ,k,K 

Moreover, we have that 



\KUcos9)\<^f-e- n . (4.3) 



H^Hp := H^bCcos^IUp^) < C rhk , K e~ n /p'. (4.4) 

These operators can be applied to functions in L p (§") or even distributions in 
£>'(§"), provided k decays fast enough - compact support will certainly work. As the 
result below shows, all them are bounded operators taking L p (S n ) — * L 9 (S"). 
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Theorem 4.2 ( |21L Theorem 3.7]) If k satisfies {4- ty , with k > max{2,n}, then, 
for all 1 < p < oo and 1 < q < oo, the operator K £n : i p (S") — > L 9 (S™) is bounded 
and its norm satisfies 

||K e) „|| p>9 <C„, fc , K (4a;„_ ie ra )-^-|)+, 

where C n ,k,n is a constant that depends only on n, fc, k, and where (x)+ = x for x > 
and (a;)+ = otherwise. 

We point out that more can be said when k has restrictions on its support. The 
result below follows from the spherical harmonics of degree L ~ 1/e being in the 
kernel of K £ „ when /-c(t) = near t = 0. 

Remark 4.3 If n(t) = /or |i| < 1, £/ien /or any spherical harmonic in TIl,, where 
L e = [e _1 — A,^ 1 ] ~ or less, then we have g e := K £jTl g = K £i „(g — P), and hence 

\\9e\\q < \\^B,n\\p,qE Le (9)p- 

Another important result for k supported away from t = and having fast decay 
is the one below, which follows directly from Theorems 14.11 and 14.21 To simplify 
matters, we will assume that k is also compactly supported. 

Corollary 4.4 Let k > max{2,n}. // the support of k is compact and does not 
include t = 0, then, for every fixed 7 in C, the function k(t) := |i| 7 >t(t) is also an 
even C k function that satisfies Moreover, L 7 K e n = £~ 7 K £i „. Finally, for real 

7, we have the two bounds below: 

\\^Ke,n6\\ P <C n , krK e-~<- n /p', 
where 5 is the Dirac distribution and thus L 7 K £i „<5 is the kernel for L 7 K e .„. 

4.2 Marcinkiewicz-Zygmund inequalities 

Marcinkicwicz-Zygmund (MZ) inequalities provide equivalences between norms de- 
fined through integrals and ones defined through discrete sums. For these were 
developed in [T71 [16l [21] . We will need to adapt these MZ inequalities to estimate 
certain sums. 

Let X C S ra be the set of centers; also, let q = qx, h — hx, and p = px '■= h/q 
be the separation radius, mesh norm, and mesh ratio, respectively. We will need a 
decomposition of the sphere into a finite number of non-overlapping regions. The 
Voronoi tessellation corresponding to X will serve our purpose here, although many 
other decompositions will work as well. 

Let i?£ be the Voronoi region containing X . Denote the collection of these regions 
by X = C S" I £ G X} and its partition norm by \\X\\ — max^ g x{diam(i?^)}. It 
is easy to show that the following geometric inequalities hold: 

h < \\X\\ < 2h and min/i(i? ? ) > c n q n . (4.5) 

Here c„ is a constant related to the volume of We will need these later. For a 
sequence space version of results below, see [T3l Proposition 4.1]. 
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Proposition 4.5 Fix ( <E S" and k > n + 2. Let := K e ^ n (r\ ■ (). T/ien, i/iere 

is a constant C = C n _ Kt k for which 



Itfellx - 5>(fle)l**(0l 



< C 



11*11 <£ 
(II^IIA)" 11*11 >£• 



Moreover, if C € X , then 

f \KMWiv)- E M«t)l^(OI 



< a 



(4.6) 



|*|| < e 



7l,K,k 



(e/\\X\\) k - n - 2 ||*|| > 



(4.7) 



Proof: The proof follows along the the lines of the one for [2lJ Proposition 4.1]. 
Therefore, we will only sketch it here, referring the reader to [21] for the technical 
details. 

The inequalities in both (|4.6|) and ()4.7j) involve bounding sums of contributions 
from each having the form 



D f 



< I \K £ (r,) - K^)\d(i(v)- 



Take C to be the north pole of the sphere and 9 to be the co-latitude. Divide the sphere 
into M ~ 7r/||;r|| bands, B m , in which (m - 1)tt/M < 9 < rrrn/M, m = 1, . . . ,M. 
Each can have non-trivial intersection with at most two adjacent bands, because 
diam(_R^) < ||*| ~ ir/M. Thus, if R^ C B m U B m+ i, then its lowest and highest 
co-latitudes satisfy (m — l)ir/M < 6*7 < 9^ < (m + 1)tt/M. As is shown in [21], for 
m = 2, . . . , M — 1 the sum of the D^ from all R^ C £? m U B TO +i is bounded above by 
the quantity, 



E 



D e < 



C 



Re CB m Ufl„ 



n,K,k 

Me 



l + t k 



dt. 



(4.8) 



If i?£ 9 Ci then dealing with the corresponding can be done by estimating the 
integral that bounds the contribution from the region R^ in the cap < 9 < 2tt/M, 



d s < a 



(Me)" 



tdt 



< 



°n,K,fc 



1 + t k ~ (Me 



(Me)~ 2 Me > 1 
1 Me < 1. 



(4.9) 



Now, let M = L7r/||*||J, precisely. Adding up the D^ for all £ € X yields the bound 
in (|4.6[) . which was implicit in the proof of [2TJ Proposition 4.1]. 

To get (|4.7p . we need to adjust M so that all i?£ ^ £ are contained in the bands 
B m U B m+ i, m = 2, . . . , M — 1. This is easy to do. Just take M = [(tt — <z)/||*||J- 
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Summing the bounded in 



and taking care of some double counting yields 



B n -R ( 



Me 

Cn ,K,,k 

Me 

Cn,K,fe 



< 



< 



t" 



oo t n 



(It 



cdt 
l + t k 

ie 

(Me)- 1 Me > 1 
(Me) k ~ n - 2 Me < 1, 



from which (|4.7p follows easily. 



□ 



Let / £ L 1 (§") and set / E := -ftT ej „ * /; the function / is not assumed to be 
zonal. We wish to estimate the difference Ex ■= |j|A||i — X^ex l/e(^)lA t (^)|- ^ i s 
straightforward to show that 

E X < E / \fM - fe(0\dKv) < sup ^(C)||/||i , 



where F B) Ar(C) : = Xfex Jr. \ K e,n(v ' C) - • C)|^m(^)j which is the quantity 

estimated in Proposition 14.51 Applying that proposition and Remark 14.31 we obtain 
the desired estimate below. 

Corollary 4.6 Let k satisfy with k > n + 2, and, /or / £ L 1 (§"), Zei / e = 

Ken * /• // <^ is the decomposition of S" described above, \\X\\ > £ and L £ = 



then 



IMIi-E l/=(OIM^e) 



<C n ,„ lfc (||Ar||/e) n 



«(t) = 0, |t| < 1. 
|| /||i, otherwise, 



(4.10) 

Remark 4.7 J// is zonal, i.e. /(£) = "0(^- C)j £«en i/ie ng/ii side zs indepen- 

dent of the variable £. j4Zso, £/ie strict inequality \\X\\ > e isn't absolutely necessary. 
The results still hold when \\X\\ and e are comparable. 

For the most part, we will use these results to bound the sums | Sfex a f/e(0|j 
under the assumption that ||«Y|| > e. Using Corollary 14. 61 for that case, we see that 



< 



< 



min 5e x M-Rs) 



Em^)i/,(oi 



(\\fe\\Li+C n ^ k (\\X\\/eT\\f\\ Ll ) 



From Theorem 14. 2 [ (|4.5[) . and /i = pa, with L e ~ e 1 and w ||Af|| > e. we have 
that 



< Cp n e-"|a| 



if «(t) = 0, I*| < 1, 
|| /||i, otherwise, 



(4.11) 
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If / is a zonal function, then, by Remark l4.7[ wc may use the || • |oo norm on the left 
above. 

We want to make the same kind of estimate, but for / being replaced by 5^, the 
usual Dirac delta function. Thus f E is replaced by K e {-) := K e , n * S(-) = A' £i „((-) • (). 
A nearly identical argument to the one used above, coupled with (|4.6p for 1 1 1 1 > e 
and the bound on j|A' £ ||i from Theorem 14. 2 [ results in 



< Cp n e- n \a\ 



The constants on the right above hold uniformly, so we thus have 

£oe* e ((-K) < Cp n S - n \aU 



(4.12) 



(4.13) 



The two bounds above are very similar and can be used in combination. They 
will be needed to complete the approximation part of the analysis. There is another 
bound, somewhat different from these two, that we will need in section [5} 



Lemma 4.8 // pq ~ || X\\ > e > and if k > n + 2, then 
max ^2 \ K e,n(£- C)l < Cn,K,M 



(4.14) 



Proof: In equation (|4.7|) . Proposition 14.51 again for \\X\\ > e, an argument similar 
to the ones used above gives us 



E \ K **(S ■ 01 < C' ntRtk q- n f \K e , n (v ■ C)\dp(v) + C',U k q- 



X3&C 



k-n-2 



Using J Bn _ R( \K e ^ n (i] ■ C)\dp(T]) < ||if e , n ||i < Cn,«,fe, -pq < 1 , and maximizing over 
C e X, wc obtain (|4"T4"]) . □ 



This estimate is more delicate than f|4. 1 3[) . because the term missing from the 
sum is A £j „(e ■ = K E , n {\). which turns out to be 0(e~ n ). For e/q small enough, 
the sum (|4.14p will be majorized by A e .„(l). This is needed as part of a diagonal 
dominance argument. 



4.3 Frames 

Wc now address the question of the frame decomposition mentioned previously. Our 
approach follows the one in |21j . As mentioned earlier, others are certainly possible. 
For this, we need a function a g C k (M) , which we may assume is even, with support 
in [-2,-|] U [|,2], and satisfying \a{t)\ 2 + \a(2t)\ 2 e 1 on [|,1]. Such a function 
can be easily constructed out of an orthogonal wavelet mask nig [2 §8.3]. In fact, if 
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rrio(£) G C fe+1 , then a(t) := m (7rlog 2 (|t|)) on [-2, — |] U [|, 2], and otherwise, is a 
C fc function that satisfies the appropriate criteria. Define b G C k (R) by 



1 Ifl < 1 



Using the properties of a we see that X)j=-oo |o(*/2 J ')[ 2 = b(t/2 J ) if i > 0. In the 
sum on the left, only terms with j > [}og 2 (t)\ contribute. Terms with j < [log 2 (i)J 
are identically 0. 

The quantity [log 2 (t) J is obviously important. On the E> n , the integer that corre- 
sponds to it is this: 

Jn := { Llog 2 °(A„)J n > 2. (4 ' 16) 

The integer j n helps us in defining our frame operators, which we now do. Let 
Aj := a(2ri~i n L„) and Bj := b(2~i~ 3n \- n ). Taking into account the support of a, 
we have Bj = J2j=o AjA* for n > 2 . For n = 1, a projection Po onto the constant 

function enters, and Bj = Po + J2'j=o AjA|. We will need the following approximation 
result concerning these operators. 

Proposition 4.9 ([21, Proposition 5.1]) Let k > max{n, 2}, and let b be defined 
by (4J5\ ), with a G C k (R). If f G L p {S, n ), 1 < p < oo, and if L > is an integer 
such that 2~ J ~i n < (L + A„) _1 , then 

\\f-Bjf\\ p <C bAn E L (f) p , E L (f) p :=dist LP (f,n L ). (4.17) 

Also, for 1 < p < oo or, if p = oo, for f G C(S ra ), we have limj^oo Bjf = f . 



Bernstein/Nikolskii inequalities. There are several inequalities that follow eas- 
ily using frames. We will give a Nikolskii-type inequality, which is a well-known 
inequality ( |15[ Proposition 2.1] and [21[ §3.5]), From our point of view, the most 
important inequality derived here is a Bernstein theorem for spherical polynomials 
[28l Theorem 2 (Eng. transl.)]. An independent proof is given in [TDl Proposition 
4.3]. For the convenience of the reader, short proofs for both are given below. 

Theorem 4.10 Let S E Hl- Then, for 1 < p, q < oo and for 7 > 0, we have 

(Nikolskii) 11% < C Piqin L n ^-^\\S\\ p (4.18) 
(Bernstein) \\S\\ h p < C nn L^\\S\\ p (4.19) 

Proof: Let 7 > and suppose L + X n < 2 J+ ^ n . From the definition of Bj, it is easy 
to see that Bj reproduces n^, and so BjS = S for all S £ Hl- By Theorem 14. 2 [ with 
K = b and e = 2~ J -^ - L' 1 , we see that < C p . q . n L n( ^^ )+ \\S\\ p , S G U L . 

Dependence of the constants on b and k disappears upon taking the infimum over 
these two quantities, yielding (|4.18|) . 



1G 



We now establish the Bernstein inequality. If 5" £ 11^, then so is L 7 S', and we have 



that BjLlS = LIS, provided L + A„ < 2 J+ ^. Using the expansion B, 7 = J2'L a AjA*, 



we see that 

J J 

US = Aj A*US = L J AjAjS. 

Consequently, we have that ||5|| ff 7 = ||L?% < E/=o II L7A j A*|| p , p ||5|| p . Applying 
Corollary 14. 41 with n(t) = |a(i)| 2 an d e = 2~ J ~ : '" for each j, then yields this: 

,/ 

oU+in)i\r< II 



\S\\ h; = < ^2(^^)^11511, 

S=0 



2^+3^ + 1)7 _ 2J*»7 
— 27_1 Ca.n.-vH'S'lIp < L~'Ca,n,i\\S\\ p , 

where again L ~ 2 J+: >™ . In the last inequality of the chain above, we can take the 
infimum over all a satisfying the requisite conditions. This yields (|4.19j) □ 



Distance estimates. Frames can be used to estimate the distance in L p (S n ) from 
the polynomials to a function in a smoother space. If / £ L p , let El(J) p := 
distLp(/,IlL). Because Bjf is a spherical polynomial in n 2 ./+j„+i, we have 

ElU) p < 11/ - Bjf\\ p , L + X n < 2 J+ ^+\ 

And because Bjf converges to / in all L p , 1 < p < oo and p = oo if / £ C(§"), we 
also have that 

oo 

E L (f) P <\\f-Bjf\\ p < ]T IIA^/Hp, 

j=j+i 

where the right side above may be infinite. Now, suppose that / = L^/i, h £ H^(E> n ), 

In that case, we have AjA*Llh = L^ (,9 ~ 7) A i A*L^/i. From this and Corollary |4~4I 
with p «-> q, we arrive at 

IIA^L^I^ = ||L-^)A,-A*L^|| P < 2-^-"<W)+>C™»)c n>fc Jft|| fl2 

Insert this in the equation above, sum the appropriate geometric series, and take 
L ~ 2 J+ ^ to get 

E 2J+jn (LZh) p < ^_ 7!a!fc! „2-^- 7 -"^-^+)(^")||/ l ||^ , 

which was essentially obtained by Kamzolov [§]. Now, since the left side above 
is unchanged if we replace L 7 by L 7 — S, S £ n 2 j+j„, we can replace by 

E 2 j+j n (\-Ph) q . Collecting these results yields the proposition below. 

Proposition 4.11 Let 7 > 0, and (3 > 7 + n(| — ~)+), where 1 < p,q < 00. If 
h £ if?, £/ien £/iere is a constant C = C n- /3 rtM such that 

E 2 j +in {V n h) p < \\(I-Bj)h\\ H?/ < C^^-^-^-^^+^E^^h),. 
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Green's functions and their properties. Let (3 > n/p'. Recall that the Green's 
function solution to \-^Gp = 5 is a kernel with an expansion in spherical harmonics 
having coefficients Gp(£,m) = (£ + A) - ^. Properties of Green's functions (pseudo- 
differential operator kernels, really) on manifolds have been studied extensively (cf. 
[8]). Our aim here is to use frames to obtain properties and various distance estimates 
that we need here quickly, and in a self contained way, for SBFs of the form (f>p = 
Gp + Gj3 * ip, where ip G L . Because the </>,g's are not in any of the Bessel-Sobolev 
spaces Hp, they have to be treated separately from the class in Proposition ^. lll above 
We will begin with Green's functions themselves. Note that kjk*Gp = L^^AjA*^. 
Since AjA* = |a| 2 (2~-'~ : '" L„), where both a and, of course, \a\ 2 , have compact support 
that excludes t = 0. we may apply Corollary 14.41 with Ej := 2~(- 7+: '"'. 

llA^G/jllp < C n ^ a e r n,v ' = C n , /3>a 2-0 s -»/f')(j+ J -»). ( 4 . 20 ) 

Thus, for (3 > n/p', the terms in ^2J^ AjA*Gp are bounded by a geometric series, 
and so the Weierstrass M test implies that the series converges in LP ' . That is, we 
have shown that when (3 > n/p' the limit limj^oo BjGp is in LP . A simple duality 
argument then shows that the kernel Gp = limj^oo BjGp in L p (E> n ). Summing the 
geometric series in KM yields \\Gp - BjGp\\ p < G2"(' 3 - ti /p')( j +^). 

These results also give us error bounds in H^(S n ). If 7 > 0, then L^Gp = G^_ 7 

and L J BjGp = BjGp r This and the estimate above imply that if in addition 

(3 > 7 + n/p', then 

\\G P - BjG \\ H , = ||G^ 7 - B,jG ^\\ p < C2-^~ n ^ J+ ^. (4.21) 

Perturbations of Gp can be dealt with, too. Let ip be in L 1 . By Theorem 14.21 
(|4.20|) and Remark 14. 3) we have that, for all j > J, 

||A,A*G /3 *^|| P < WA^GpW^Eu+jnWi < C2-^-^ + ^E 2 ., +]n ^) l 

Summing a geometric series and using (|4.2ip , we arrive at the following bound. 

Proposition 4.12 Let 7 > ; (3 > j + n/p', ej = and let ip € L 1 be a zonal 

function. Lf (ftp = Gp + Gp * ip, then <pp <G and there is a constant C = G„^ j7ia; 
which depends only on n, (3, 7, and the function a, such that 

E 2J+jn (U<f>p) p < \\(L-B, ; )cf,p\\ H , < C^^l + e^'E^i^y^-^'. 

(4.22) 

4.4 Approximation analysis 

The task at hand is to estimate the norms ||(/— Bj)g|j^fp/|a| p , where g G Gx,4>- Our 
approach will be to carry this out for p = 1 and p — 00, then use the Ricsz-Thorin 
theorem to obtain the result for all intermediate values of p. 



18 



The easier of the two cases is p = 1. Since g <E Gx,<j>, we have g = X^ex a ?<M(') '£)• 
Again, let = 2 From the triangle inequality, the rotational invariance of 

the norms involved, and Proposition 14. 1 II and Proposition 14. 1 21 it follows that 

\\(I-Bj)g\\ H i < {aUWil-Bj^Wni 
The p = oo case requires using frames. Again, we have that 

oo 

\\(I-Bj)g\\ HT < J2 IIA^L^I^, 

j=J+i 

where AjA^g = £ £e jc a e A,-A]L^((-) ■ £). By equation gH|), with / = L^, A £j ,„ 
corresponding to n(t) — \a(t)\ 2 , h> sj > Sj, all j > J, and L e ~ 2 J +- ?n , we have 

llAiAjLZflHoo 



<Cp" £ T"l«|o0^n(L^)l 



IIA.A^L^IU < Claire* 



nJ—7—n 



By Proposition 14.111 and Proposition 14.121 with J there replaced by j, p — oo, we 
have 

E 23+3n (L^)i e , 
(l + e J n A 2J+J „(^) 1 ) = G^ + G^*^. 

Since E 2 j+j n (/)i < E 2 j+j n (/)i when j > J, in the inequality above we may replace 
the distances with respect to 2 J+J ™ with ones with respect to 2 J+ -'" . Doing so and 
again summing a geometric series, we obtain 

Applying the Ricsz-Thorin theorem in conjunction with the bounds above, we com- 
plete the approximation part of the problem: 

Theorem 4.13 Let 7 > 0, 1 < p < 00, (3 > j + n/p', £j = 2"( J+J "). If h x > £j and 
if 9 £ Gx.4,, then 

\\{I- Bj)g\\ H P , «_ 7 _„ /p ' f S 2 .7 +J „ (L^)i e ^ , 

|o| p - P J I (l + ^Wi) <f> = Gp + Gp*i). 

(4.23) 



5 Stability 

The problem that we address here is estimating the norm \a\ p in terms of the L P (S") 
norm of g, where g(x) = X^e x a £0( x ■ an d <fi E L p is an SBF. Specifically, we wish 
to estimate the p-norm stability ratio 

Hp 

r g p := max — — 

03S*O \\g\\ p 
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which we defined in This quantity exists and is finite because the set {</>(x • 

0}?ex is a linearly independent, finite set of functions. The quantity xq^ v provides 
a measure of the linear independence of the set, albeit one that scales with the norm 
of cj). Once <f> is fixed, it depends completely on the geometry of X. 

For a continuous SBF cf>, this is related to the stability of the interpolation matrix 
for cj> and X. However, we are only assuming that <f> is in L p , and thus evaluating <f> 
on X is meaningless. Even so, using a smoothed version of <fi allows us to connect 
the two concepts. 

5.1 Stability ratios and interpolation matrices 

Let k > be in C fc (R), k > n + 2, and let it satisfy (|4.ip . Of course, since n is 
not identically 0, we also have that there is some open interval on which k > 0. 
Consider the corresponding operator K e .„ = /t(eL n ) and its kernel K e n . To smooth 
<?(x) = X^£<ey a £0( x ' £)j a PPly K e .n to both sides. Doing this yields 

g e (x) = K s , n #(x) = a ? K e,«0(x • (5-1) 



We want to relate rg jP to quantities in a standard SBF interpolation problem on 
X involving <f> e . The function (f> E is a spherical harmonic, with nonnegative Fourier- 
Lcgcndrc coefficients, whose degree depends on the support of k. It is thus a positive 
definite function on S n , but not an SBF. 

The interpolation matrix corresponding to (f> e is 

A e = y> E (r) ■ Oh,nex- 

Later, as a by-product of our analysis, we will establish the invertibility of A £ , pro- 
vided e satisfies certain conditions. When e is sufficiently small, one can also establish 
it by using a result of Ron and Sun [571 Theorem 6.4]: Let X C §" be fixed and let if) 
be a positive definite function, but not necessarily an SBF (i.e., some of coefficients 
ip(l) may vanish). Then, there is an integer jx,n such that the interpolation matrix 

will be positive definite if the set of integers on which ifj(£) > contains at least 
jx,n consecutive even integers and jx,n consecutive odd integers. With our assump- 
tions on k — in particular, that k is not identically - it is clear that for sufficiently 
small e there are arbitrarily large sets of consecutive integers for which (f> £ (£) > 0. 
Thus A e is (strictly) positive definite, and hence invertible, for all such e. 

Our approach will again be to use the Riesz-Thorin theorem. Let y e := g e \x, the 
restriction of g £ to X. Using (|5.1[) . we can interpolate g e on X: 

lJe = A E a , A E = [<p £ (-q ■ Oh,nex, 
Solving and taking the i 1 norm, we see that 

Hi^il^Hxi^u, iifeii = X>=(0I- 
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By our assumptions on k and by (|4. 1 1[) . we have that \y e \i < C n ,n,kP n £ "HfflU 1 - 
Consequently, for </> £ L 1 we have that 

rg,i<C K>n!k p n £- n \\Aj%. 

Similarly working with p = oo we obtain 

Moo < Pr^loolyeU > \Ve\oo = max{\g £ (£)\} < ||g||oo- 

Recall that A7 1 is a self-adjoint matrix, and that for such matrices the p = 1 and 
p = 00 norms are equal: HA" 1 ^ = H-A" 1 !!!. Hence, for (f) £ C (p = 00), we obtain 

Applying the Riesz-Thorin theorem to these bounds yields the following: 
Proposition 5.1 Let e < \\X\\ and let <j) £ LP . Then, 

^,v<C l Jl kP n/p e-^\\A-%. 

5.2 £ l stability estimates for interpolation matrices 

The estimates we need next are for 1 1 ^4 ~ 1 1 [ 1 , and the approach we take to get them 
will depend on cf> and the behavior of the 4>(£) 's. We will first deal with the Green's 
function case, in which 4>(£) decays algebraically After that, we will deal with the 
case in which </> is C°°, and (f>(£) has very fast decay. 

5.2.1 SBFs that are perturbations of Green's functions 

A straightforward way to estimate the 1-norm of the inverse of a matrix is to use 
diagonal dominance techniques, if the matrix is amenable to them. To that end, split 
an n x n matrix A into its diagonal D and off-diagonal F, so A = D + F. We then 
have the following standard norm estimate, whose proof we omit. 

Lemma 5.2 Lf D is invertible and ||-D _1 -F||i < 1, then A is invertible and < 

WD-^ui-WD-^FhyK 

We can apply this to A e . The diagonal part is D = <fi £ (l)L, and so j|_D _1 ||i = 
e (l) _1 and = <t^^ ( IL ) 1 1 1 | [ 1 - Since the 1-norm of a matrix is the maximum 

of the 1-norms of its columns, our condition becomes 

Miy'WFh = M^r 1 ™? E \Mv-o\<i- (5-2) 

We now want to deal with a special <p £ , which is not necessarily generated by an 
SBF (p. Let ip be a zonal function in L 1 , so that 

1=0 XnUJn 
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We will assume that 1 + t/j(£) > for all £ > and that k has support in \t\ £ [1, oo). 
Take <fi e = K e>n + K £>n * ip, where K s>n is the kernel for the operator k(eL). In 
addition, define ip E = K e ^ n * -0. Since 4» e {£) = n(e(£ + A„))(l + 4>{£)) > 0, we see that 
<fi e is a positive definite spherical function, but not an SBF. Using (14.14j) yields 

E \Mv-0\ < E \Ke,n(vt)\+ E 

XB&ri X3^ v X3^ v 

icx 

Thus, from this and equation (|4.1ip . with n(t) = 0, \t\ < 1, we have shown that 

E \Mvt)\ < C n ,K, fe (g-" + p n £- n i?L e Wi), L s = Ll/e-AnJ. (5.3) 

Thus we have bounded the sum involved in the diagonal dominace condition (|5.2|) . 
Next, we will deal with e (l). We have the following chain of inequalities: 

e (l) - K s>n (l) + K E , n * if>(l) 

OO 

= E K ( £ ( £+A «))(i+V'WK 

oo 

> c E + A »)K = 

where Co = min^>o(l + V'M) > 0. (This is true because tp € L 1 implies that ijj(£) — ► 
as £ — * oo.) Furthermore, it is easy to see that 

°o />oo 

K e , n {l) = E + X n))d" ~ e~ n / n(t)t n - l dt . 

t=o w , 

>0 

Thus, </> e (l) > C'' fc £~™. From this and (|5.3p . we arrive at the bound below: 

p-^lli < C n , K , k ((e/q) n +p n E Ls {iP) 1 ), L e =[l/e-X n \. (5.4) 

By choosing e < q sufficiently small, we can make C n ,K,,kP n EL r (ip)i less than 1/4, 
since El s (iP)i — ► as L £ — > oo. At this point, the choice of e depends only on ip and 
the mesh ratio p. If necessary, we may then choose e smaller still in order to force 
the first term on the right to be less than 1/4. With this choice of e, which depends 
on p, n, k and k, we obtain 1 1 ^ZH> 1 ^f 7 ^ 1 1 1 < 1/2. By Lemma 15.21 we get the bound on 
|| A" 1 ||i below. 

Proposition 5.3 Suppose that k has support in \t\ e [1, oo). Let <f> £ = K e ^ n +K etn *^, 
where ip £ L 1 is a zonal function satisfying 1 + ip(£) > for £ > 0. Then there are 
constants c and C , which depend on ip, on p, n, k and k, such that whenever e < cq 
we have {{Aj 1 ^ < Ce n . 
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The proof above required conditions on the support of K in order to deal with the 
perturbation generated by ip. If ip is 0, then there is no need for such restrictions. 
Also, the term involving p is gone, and it is no longer involved in determining c and 
C. We collect these observations below. 

Remark 5.4 If tp = ; then Provostion \5.3\ holds without restriction on the support 
of K, and neither c nor C depend on p. 

We now take an SBF <j> of the form <fi = Gp + Gp * ip, where Gp is the Green's 
function for YJ 3 and ip € L l . Our aim is to establish a bound on the stability ratio 
for such <f). 

Theorem 5.5 Consider the SBF cp = Gp + Gp*ijj, where Gp is the Green's function 
for and ip E L 1 . Let X be a set of centers with separation radius q and mesh 
ratio p. Let Q — Q^.x be the corresponding SBF network. Then there is a constant 
C = G{n, </>,/?) such that the stability ratio of Q satisfies 

rg, P <Cp n/v q n/p '- r3 (5.5) 

Proof: Since we are assuming that <f> is an SBF, the coefficients of the L 1 function 
tfj must satisfy 1 + > for all I > 0. Assume k satisfies (|4.1j) and has support 
in \t\ £ [l,oo). The corresponding <f> e is just 4> e = K e „0 = K £t7l (Gp + Gp * ip). 
By Corollary 14.41 we have that K e ^ n Gp = e^K £) „ = «(eL), where k{t) = |t| _/3 K(f) 
satisfies (|4.1[) . From this, we have that (j) £ = e^(j) e . If we let A e be the interpolation 
matrix for e , we see that A e = e@ A e . The function <fi e satisfies the conditions on the 
corresponding function in Proposition [5T3I Thus, by choosing e < cq, we have 

|| 1 1| ! = e-^H^r 1 ||i < Ce n ~P. 
From Proposition we obtain 

r e , P < C^p^e-^WAjX < C'p" /p e" /p '- /3 , 
Choosing e as large as possible, namely e = cq, we have 

r e , P <Cp"/ p 9 "/ p '- /3 , 

where the constant C = C(n, k, k, <fi,p, ft). By taking the infimum over all K, p and 
k, we reduce the dependency of C to C — C(n, <f>, 0). This completes the proof. □ 

5.2.2 Infinitely differentiable SBFs 

Let </> be infinitely differentiable SBF. The fast decay of the Fouricr-Lcgendre coeffi- 
cient 4>(£) requires a different approach to bounding xg than the one used to obtain 
Theorem l5.5l As before, we let A e be the N x TV interpolation matrix for (j> £ = K £: „</j. 
In addition, we will let A be the corresponding matrix for </>. By standard matrix 
estimates, the norm 1 1 v4.~ 1 1 1 i satisfies 

\\A-% < N^ 2 \\A-%. 
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Since A e is a positive definite selfadjoint matrix, the norm HA" 1 ^ is equal to the 
reciprocal of A m i n (^4 e ), the smallest eigenvalue of A E ; that is, H^" 1 ^ = 1/A m i n (^4 e ). 
We will begin by estimating this eigenvalue. In preparation for this, we define the 
quantity 

min (L) := min > 0. (5.6) 

a<e<L 

where the strict positivity follows from <f> being an SBF. 

Proposition 5.6 Let re > be in C k (R), k > n + 2, and let it satisfy In 
addition, suppose that supp(re) C [—2,2] and that re < 1. Then, there are constants 
c = c n _ Kt k > and C = C„. K ,fc > such that for all e < cq, 

A m in(A) > X min (A £ ) > C4> min (L e/2 )£- n , L e/2 := [2/e - A„J. 

Proof: Using the Raylcigh-Ritz principle, we thus have 

P^Ha 1 = \mn(A e ) = min a* A £ a. 

aeC N 

where A £ = [4> £ (v ' 0]$,riEX ■ Because <p e is a (positive definite) zonal function, we can 
use its expansion in spherical harmonics to represent A m i n (A e ) via 



A mi „(A £ ) = min I V V K ((i + \ n )£)m 



e=o m=i 



2 



(5.7) 



Since the support of re is [—2,2], the sum above cuts off at L e / 2 := [2/e — A„J. 
Consequently, we can bound below X m i n (A £ ) this way: 



e/2 



Amin 

(A E ) > (j} min (L £ / 2 ) rain > > n((£ + X n )e) 
a€C N 1 ' — ' L — ' 



g=0 m=l 



X m in{[K e ,n{i ■ V)]) 

Note that \ min {[K e , n {i ■ rj)]) = \\ [K e>n (£ ■ vT 1 ^ 1 < II [K e ,n(Z ' »?)]~ 1 llr 1 > because 
II-BII2 < ll-^lli f° r ai l selfadjoint £?. The existence of c and C and their dependencies, 
along with || [K e , n (£ ■ ^)] _1 ||i < Ce" for e < cq, follow from Proposition 15.31 and 
Remark 15.41 Finally, applying the Rayleigh-Ritz principle, (|5.7p . and < re < 1, we 
have that \ m [ n (A) > X m i n (A e ). This finishes the proof. □ 



There are two immediate consequences that follow from Proposition l5.6l The first 
is a bound on the stability ratio in this case. 

Theorem 5.7 Consider the SBF <f>, where <j> is assumend to be infinitely differen- 
tiable, and let X be a set of centers with separation radius q and mesh ratio p. Let 
Q = G& x be the corresponding SBF network. Then there are positive constants 
C = C n _ K k and c = c„. K ^ such that the stability ratio of Q satisfies 

n(l/p'-l/2) 

rg, P < Cp n/P j-— -, where L cq/2 = \2/{cq) - A„J 



24 



Proof: Since 1 1 1 1 1 1 < iV 1 / 2 1| Aj" 1 1| 2 , Proposition I5H implies that for e < cq, 

N X/2 £ n 



\A- l h<c n 



By Proposition [O] we then have that 

pn/p e n/p' 



tn(i £ /2) 



Noting that N ~ q " and choosing e = eg, which is as large as possible, we obtain 
the desired inequality. □ 

The second consequence is a new stability estimate for interpolation via a C°° SBF 
<j>. Again, let A be the interpolation matrix for <j> on the set X. By Proposition 15.61 
!l^4 _1 ||2 = A m in(A) _1 < Ce n /(j) m i n (L e /2)- Taking e = eg, we obtain a new bound on 
the norm of A~ 1m . 

n 

\\A-%<C 1 -^ r -. (5.8) 

<Pmin{J-> cq / 2 ) 

6 Bernstein inequalities and inverse theorems 

In this section, we will discuss both direct and inverse theorems for approximation 
by SBFs. For an overview of these notions, see [5]. 



6.1 Bernstein inequalities 

Bernstein inequalities are a primary tool in obtaining inverse theorems. In the intro- 
duction, we gave a strategy for obtaining Bernstein theorems. We have completed 
the preparation required to state and prove them. Our first result is for SBFs that 
are perturbations of Green's functions. 

Theorem 6.1 Consider the SBF cp = Gp + Gp * ip, where Gp is the Green's function 
for L 73 and tp € L 1 . Let X be a set of centers with separation radius q and mesh ratio p, 
and let Q = Q^.x be the corresponding SBF network. If 1 < p < oo, < 7 < fl — n/p' 
and g S Q , then 

IMk < Cq-*<\\g\\ p . (6.1) 

Proof: Recall that ||g||#-p < ||-Bj<?|| ff »> + — B>j)g\\ H v, where Bj is the frame 
reconstruction operator defined in section Of course, from (|4.17p . this operator is 
bounded independently of J From the polynomial version of the Bernstein inequality 
in P~T§)) , we have that \\Bjg\\ H P < C2~< J \\ Bjg\\ p < C2> 7 ||. 9 ||p, which implies (fT^j) . 
Inserting the approximation estimate (|4.23|) and the stability-ratio estimate (|5.5p into 
(fTT2|) yields 

\\g\\ H * < (CT< J + C'2-^- n ^ J q^'-(\l + E 2 . J+]n ^) 1 ) \\g\\ p 
< q-i (C{2 j qy + C'(2- J q)^- n ^(l + |M|i)) \\g\\ P 
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The integer J is still a free parameter. Choose it to be J = — log 2 (q). The Bernstein 
inequality (|6.1[) then follows on noting that q < ir, (3 — 7 — n/p' > 0, and ||V'||i is 
finite and fixed. □ 

Up to a point, an SBF <f> <G C°° is handled in the same way as one related to a 
Green's function. In particular, using the argument above, coupled with the approx- 
imation estimate (|4.23[) , with (3 = 7 + n, and the stability estimate in Theorem 15.71 
we obtain 

\\9\\ m <CL> (l + CV"( g L)"(W-V2) L ~*~* ) ? T M * )l ) \\g\\ p , L = 2^», 

(6.2) 

where L cq / 2 = [2/cq — A„J . Because (f> € C°°, it is in for all [3. The inequality 
thus holds for all (3 > j + n/p'. The object here is to find a constant L = aq^ 1 , where 
a is independent of q, such that the ratio on the right above is bounded. The other 
terms will be controlled easily in that case. To obtain a simple, applicable condition, 
we need the following lemma. 

Lemma 6.2 Let < j-i{£) < a(£) be eventually decreasing sequences. Assume that 
for every a > there is an integer mi = mi(a, a) > such that £ a o~(£) < o~(2~ mi £). 
If in addition for all £ sufficiently large there is an integer m2(tt, £t, cr) > such that 
a(2 ni2 £) < Cfj, t<7 fj,(£), then with m = mi + m2, 



1 



J2 £ a o{£)<C^2- m L- 



=2 m L 

Proof: Let mi = mi (a + 2,4>). then 



00 



e=L £=l e=L 

Replace L by 2 m L in the inequality above, so that the sum on the left above is 
bounded by {2 m L)- 1 a{2 m2 L) < C ' Ml0 -2- m L' 1 fi{L) . Dividing by fx(L) yields the de- 
sired inequality. □ 



Lemma 6.3 If there are two sequences /i(£) and o~(£) that satisfy the conditions of 
Lemma 1 6. S\ and in addition satisfy fi(£) < cf>(£) < <r{£), then there is an integer 
m = m(/3, (f>, n) such that for all L sufficiently large 

< C M , n 2-™ L -\ (6.3) 

(Pmin(L) 

Proof: Because is a C°° SBF, the error i?i(L^</))i satisfies 

El (l^)i < u n E L {i^ )oo < £ (£+A - ) ^ (A " )(1) 0ffl < ^ E ^ +n - x m 
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where we have estimated factors independent of 4> to get the term on the right. 
Applying Lemma ROl then completes the proof. □ 



Putting all these results together leads to this theorem. 

Theorem 6.4 Let <j> be a C°° SBF. If there are two sequences fj,(£) and a(£) that 
satisfy the conditions of Lemma \6.2\ and in addition satisfy fj,(£) < <$(£) < o~(£), then 
for every 7 > Bernstein's inequality, 

WgWw; < C0 :7 ,p<7 _7 ||.9||p, 

holds for all g £ G$,x, 1 < p < oo. In particular, it holds for the Gaussians, 
multiquadrics, ultraspherical generating funstions, and the Poisson kernel. 

Proof: To get the inequality itself, use Lemma I6T41 with /3 = 7 + n > 7 + n/p' . The 
statement concerning the list of functions may be established by checking that upper 
and lower bounds given in section [3] for each function satisfy the conditions on /j,(£) 
and a{£). □ 



6.2 Direct theorems 

In [T5l §4], we used a linear process to estimate the distance distLp(/, G<p,x), given 
that <f> is a continuous SBF and / € L p . In several important cases, including the 
Gaussian, the process produced a near-best approximant. We will use a similar 
process here for an SFB of the form <fip = Gp + Gp *tp, tp £ L 1 , again obtaining the 
corresponding distance estimates. Such SBFs are at least in L 1 , but they might not 
be continuous. Our approach also makes use of recently developed positive- weight 
quadrature formulas for S n , introduced in [T7] and further developed in |JT]. We 
remark that a version of Theorem 16.81 with the conditions on <f> given in terms of 
sequence spaces involving the </>(^)'s, was established in [T3J Theorem 3.1]. 

The general framework is this. Let (j> be an SBF, so that the Fourier-Legcndrc 
coefficients cj>(£) are positive for all £. Define </> -1 to be the formal expansion 

An Ulm 

£=0 

This expansion will converge in a distributional sense if the </>(^) -1 grow polynomially 
fast. Otherwise, i.e. for faster growth, the expansion is purely formal. Since we are 
using it in connection with polynomials of finite degree, this is not a problem. 

For every spherical polynomial S G Hl, we can use </> -1 to define an inverse for 
the convolution operator S — » (/> * S € ITx; namely, the expression * S, which is 
defined by the expansion 
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which is just the convolution of S with the polynomial Yli=o a + ^" l Pt ■ 

Suppose that 5 is a spherical polynomial for which deg S + A„ < 2 J+: > n . By 
Thcorcm l4.10[ we have that BjS = S. In addition, S = <f>* (j)^ 1 * S. Combining these 
two then yields 

S(x) = B J (f>*(f>- 1 *S = f {Bj<f>){x- v ){<f)- 1 *S){ji)dfi{?i). 

The kernel Bj4>(x-t]) is a zonal polynomial with degree less than 2 J+: > n+1 . In addition,, 
(j)^ 1 * S is a spherical polynomial of degree 2 J+J ' n_1 . Thus, the integrand above is a 
polynomial of degree less than 2 J+J " +1 + 2 J+jn ~ 1 < 2 J+j "+ 2 . 

We will discretize this integral by applying the quadrature formula in J5TJ §4.2]. 
Let X be a set of centers, with q, h, p, and X being the separation radius, mesh 
norm, mesh ratio, and Voronoi (or similar) decomposition, respectively. Take L > 
be an integer. There are positive weights c*, £ <E X and a constant s„ > (cf. [2U 
§4.1]) such that 

/ f(r])d»( V ) = (6-4) 

holds exactly for polynomials in 11^, provided that /i < jS~ 1 (i + A ri ) _1 . The weights 
behave like Qr = O (h n ), where the constants hidden by "big" O are dependent only 
on the dimension n. Applying the quadrature formula to the integral representing S 
yields 

Of course we are assuming that h ~ 2~ J . Let Q : IIl — ► be given via 
Q e S(z) := X! c €^-0(^ _1 *5)(0, 

and let g = QgS, where Q is used because of the operator's relationship with quadra- 
ture. The difference between g and S is thus 

g-S=J2 <*(I - BjM(-) ' CXr 1 * = (I- Bj)g. 

We now want to estimate the norm of the difference g — S = (I — Bj)g in 
terms of || c^^ 1 * S'llp- It is important to note that the norm * <S|| P depends on 
the degree of S and on <j>. We will deal with it later. 

The easiest way to estimate \\g — S\\jjp is to employ Theorem l4.13l where the norm 
ratios ||(/— Bj)g||jyp/|a| p have been estimated. Thus, the task to be accomplished is 
to relate |a| p to ||</> -1 * S\\ p . To do this, we will again use the Riesz-Thorin theorem. 

First of all, we have that a^, which is the coefficient of </>((•) • £) in g, is given by 
a c = c^cf)- 1 * S)(£). Thus, |a|oo = max^x^K^ 1 * S)(£)\. Since c 6 = 0{h n ), the 
bound |a|oo < Ch 71 ^ 1 * S\\oo holds. 
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The p = 1 case requires more work. Now, |a|i = X^ex c d(^ 1 * S)(£)|. Since 
c ( = O (h n ) < C nP n q n < C£p n min 5eA - n(R t ) < C'^p n p(R{), we have 

Mi < c>" J2 aWIW-" 1 * s)(0\ < —f-U~ l * • 

The right side above follows on applying the polynomial version of the Marcinkiewicz- 
Zygmund inequality from [211 Theorem 4.2], with S = 1/4, to bound the sum in the 
middle by (5/4)||0" 1 * S\\i. The Riesz-Thorin theorem then implies 

\a] p <C n #p n / p h n ' p '\\<t>- 1 *S\\ p . 

Combining this with the estimate (|4. 23[) . where h ~ e.j = 2~( J+: '"' and noting that 
g — S = (Qg — I)S, we obtain the following result. 

Lemma 6.5 Let 7 > 0, 1 < p < 00, (3 > 7 + n/p', h ~ 2~^ J+ ^ n ' . If S is a spherical 
polynomial of degree 2 J+ 3™~~ 1 or less , then 

kq. - m« < v^ir' . n, { * =%' + e„*. 

The <f)fj case. We will now focus on the <^a's. Our immediate concern is esti- 
mating \\4ip 1 * Slip. 

Lemma 6.6 Let 1 < p < 00, [3 > 0, ip € L 1 , and S <G U L . If <f>p = Gp + Gp * V, 
i/ien i/iere is a constant C = C n , P ,^, which is independent of (3, L, and S, such that 
this holds: 

II^SH^a^llSII^. (6.5) 

Proof: Note that <^>7 * S — (L^ipp)^ 1 * L^S. The kernel Gp is a Green's function 
for and so \Jfyp = 5 + 6*tp = d + ip, which is to be regarded as a distributional 
kernel. Finding (L^)" 1 !.^ requires solving Lfyp * T = T + ip * T = L%S for T in 
Hl, which can be done directly, coefficient by coefficient. The solution T is of course 
unique. 

There is another way to look at this equation, in an LP setting. Suppose that 
we want to solve Hf := f + tp * f = h in LP, for 1 < p < 00 and in C (for 
p = 00). The operator norm for / — * rj) * f is H^Hi- By Theorem 14. 91 we have that 
||f/> — B,/-0j|i — > as J — > 00. It follows that the convolution operator with kernel ip is 
the norm limit of finite rank operators with convolution kernels, Bjip. The operator 
ijj* is therefore compact on all LP and C; hence, Hf = / + ip * f has closed range on 
these spaces. Moreover, a simple coefficient argument shows that ker(if) = {0}. The 
Fredholm Alternative [H § VII. 11] then implies that ker(-fP) = {0}, so H~ x exists 
and is bounded on all LP and C. Since <j>^ * S = iT _1 L^S, we have that 

llV^Silp^lll^UISII^. (6.6) 
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We emphasize that \\H 1 || p is independent of /?, L, and S. It depends only on p, n, 
and V'- Consequently, C U)P ^ = and (|6.5[) holds. □ 

These lemmas lead to the following two direct theorems, the first for S € Hl and 
the second for / £ H^. 

Theorem 6.7 Let 1 < p < oo, 7 > 0, and [3 > j+n/p' . If S is a spherical polynomial 
of degree 2' /+Jn ~ 1 or less and if h = pq ~ 2~' / ~-'™, £/ien we ftai/e /or 4> = <f>p, 

h^\\S\\ H$ , (6.7) 

Proof: The two lemmas, when applied to <f)p, yield 

II (Qe - I)S\\m ^ C n , 0tl , M p n h^\\S\\ H *. (6.8) 

The result follows on observing that dist#p (S, G<f>p,x) < ||(Qe — I)S\\h p ' Note that 
the dependence of C on the particular frame operator disappears on minimizing the 
constants involved over all functions a. □ 



Theorem 6.8 Let 1 < p < oo, 7 > 0, and /3 > 7 + n/p'. If f e H p , then for 
4>fi = Gp + Gf3*ip, V G L 1 , 

dist H p(f,g^ x ) < C f3n ^ p ^h f} -~ t p n \\f\\ H P . 

Proof: Let 2~ J ~' Jn ~ h and choose S to be the polynomial S — Bjf; note that 
QgS £ G<j>a,x- From these choices and (|6.8[) . it follows that 

IIZ-Qc-SHfl* < ||/-B J /|| fi3 P + ||(Q { ,-J)S|| fl * 

< 11/ - BjIWh* + h^p n C , n JBjf\\ H *. 

By Proposition 14. 1 fl with p — q, we have and 

11/ - Bj.f\\ m < Cp, w 2-^-^ J +^E 2 ., +]n {Lif) p < C^, n , a h^\\f\\ H P. 

From Proposition 14. 9[ we easily see that ||Bj/|| ff p < C J g j -y jTlj0 ||/||.H-j>. Combining all 
of these inequalities establishes that 

11/ - QgSWn* < C^ w4 ,p n h' 3 -''\\f\\ H , (6.9) 

Since distjyt>(/, G</>,x) < 11/ — QgS\\h*> an d since the distance itself doesn't depend 
on the particular frame function, minimizing over the a yields the result, with the 
constant independent of a. □ 



30 



The C°° case. The case in which the SBF <f> is C°° was in large part done 
in [15] . However, some adjustments need to be made because the estimates in that 
paper did not involve H?. One difference is in estimating the norm \\<f>~ 1 * S\\ p . 

Lemma 6.9 Let 1 < p < oo, S > 0, L > an integer, and S £ Hl- If <t> £ Hf is an 
SBF, then there is a constant C = C n , depending only on n, such that this holds: 

T n \h~-\ 

U- l *S\\ p <C„ — \\SWm, (6.10) 

where L^ rain (L) = mxn <i< L {l + X n ) s 4>(£). 

Proof: We begin by estimating *S , || P . The case in which <fi 6 Hg was essentially 
done in the proof of [15j Theorem 4.1]; the result, which makes use of the Nikolskii 
inequality (|4.18p . is the following. If S € IIl, then the Nikolskii inequality implies 
that 

H0- 1 * S\\ p = ll(L^)" 1 * L s n S\\ p < C n L»<W>+ ||(L^)-* * L s n S\\ 2 ■ 

At this point, we simply use the 2-norm estimate done in [151 Theorem 4.1] and a 
second application of (|4.18|) to get 

ll(L^)" 1 * L s n S\\ 2 < nKiniLT^KSh < C n L n ( l *-i^([y> m jL))-^LiS^ . 
Putting the two inequalities together completes the proof. □ 

Let (j> € C°°. We can now estimate the distance of 5* € II ^ to G<f>,x, m terms 
of ||S||ijp, where <5 > 7 + n/p'. In Lemma r6.5[ let f3 = 5 + n/2. Apply Lemma T6.91 

noting that L < 2 J +^- 1 < /i" 1 implies < L n/2 < h~ n / 2 to get this: 

dist H ,(S,G^ x ) < ||(Qs - I)S\\X! ^ Cn, pP n h 5 -~< 2J+ JZ} " ^\\S\\ H * . (6.11) 

Theorem 6.10 Let 1 < p < 00, 7 > 0, S > 7 + n/p' , and <fi £ C°° . If there is an 
integer m = m{8 1 <p) > such that 

£2"^(Ln + " /2 ^)l 

sup = <C„ i; „.5 (6-12) 

holds, and if S e H L , with L < 2 J+ '^~ 1 ~ m and h ~ 2~ J -^, then 

distf^S", Q<t>,x) < C I „,„,p^, 7 /i' 5 ~ 7 p"||5'||_f/|' ■ (6.13) 
In addition, for f S H?, we have that 

dist ff p(/,^,x) < C m ,n, P ,^h s -^p n \\f\\ H P. (6.14) 

Finally, these estimates hold for Gaussians, multiquadrics, ultrasherical generating 
functions and Poisson kernels. 
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Proof: If (16T21) holds, then, since L 5 n (j> min (L) > L^ min (2 J+ J"^ m ), it follows that 



\ — — ^ (_ 

L^ mi n(i) " l40«in(2 J+J '»- OT ) 



and (|6.13|) follows form this and (|6.11|) . One can establish the H p distance estimate 
(|6.14p below using a proof virtually identical to that for Theorem 16.81 Essentially 
the same argument used in section 16.11 can be used here to show that Gaussians, 
multiquadrics, etc. satisfy (|6.12|) . and so the estimates hold for them, too. □ 



6.3 Besov spaces. 

In this section, we review the definitions and basic facts regarding Besov spaces on 
These spaces, which will interpolate between L p (E> n ) and H p , are defined in [34] . 
Other, equivalent definitions of Besov spaces on §™ are given [20] . Below, we will 
make use of a general construction general construction found in [5] Chapters 6] to 
characterize these spaces in terms of spaces of SBF networks, Gc/>,x- 

There arc two ingredients. First, wc need to introduce certain sequence spaces. If 
r > and < r < oo, wc define for a sequence a = {a n }^Lo of real numbers, 

sup2™>„|, 

^ n>0 

The space of sequences a for which ||a|| T:r < oo will be denoted by b T: ,.. 

The other ingredient in the definition of Besov spaces is a /Sf-functional [3[ Chap- 
ter 6]. For S, 7 > 0, 1 < p < oo and / £ L p , the if-functional for L p and H p is given 

by 

ic 7 (pj,6):= inf{||/ -g\\p + ^(\\g\\p + \\g\\^)}- (e.ie) 

gent; 

If r > 0, < t < oo, r < 7, wc define the class of all / £ L p for which 

||/||r. 7 ,r, P := ll/ll P + l|{iC 7 (p,/,2- n )}~ || T , r <oo. (6.17) 

to be the Besov space B r T As we shall see, other than the requirement r < 7, the 7 
dependence will disappear from the characterization of the space, so it isn't necessary 
to keep it in designating the space. 

An important problem in approximation theory is to characterize Besov spaces 
using degrees of approximation of functions. We recall the results [3j Theorems 7.5.1 
and 7.9.1], as it applies in the context of the present paper. 

Proposition 6.11 Let 1 < p < 00, 7 > 0, and let {Vj}°^ , with Vo — {0}, be a 
nested sequence of finite dimensional linear subspaces of L p , p < 00 or C , p = 00 
Suppose that for j = 1, 2, • • • , one has both the Favard (Jackson) estimate 

diflt L p(/, Vj) <C2-"(||/|| P + 11/11^), (6.18) 



if < t < 00, 
if r = 00. 



(6.15) 
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for all f G HP and the Bernstein inequality 



WqWh* <C2^\\g\\ 



(6.19) 



Then for <r < 7, < r < 00. / G £?£ p i/ and onfr/ if {disttP (/, Vj)}jt G b T , r . 

Proof: This is just [H Theorem 7.5.1], with the sequence of spaces satisfying all 
requirements in listed in [21 (5.2), p. 216], except possibly density. This requirement 
is in fact satisfied if the Favard inequality (|6.18[) is satisfied. To see this, note that 
contains all of the spherical polynomials, which form a dense set in L p , 1 < p < 00 
and in C. The Favard inequality (|6.18p then implies that the UjVj is dense in 
and therefore in LP , 1 < p < 00, or in C. □ 

In the important case when Vj = n 2 j , Proposition 14.101 gives the Bernstein es- 
timate, while Proposition 14.111 provides the Favard estimate. In addition, since the 
criterion that {distj>(/, II 2 j)}^_ G b T r docs not depend upon 7, it follows that the 
Bcsov spaces B r T p are independent of the different choices of 7 > r in their definition. 
This is why we don't need to include the parameter 7 to index these spaces. 

Remark 6.12 The polynomial characterization of ' B r T is precisely the one given in 
\2(A Proposition 5.3], so that the "needlet" definition \2(A Definition 5.1] is equivalent 
to the one above. (See also ]18jj .) The needlet definition is itself known to be equivalent 
(cf. ]20jl ) to that given in \34)j - It follows that all three are equivalent. 

Using the proposition above, one can also characterize Besov spaces using a va- 
riety of spherical basis functions. To do this, we must first have an appropriate 
nested sequence of sets of centers. By Proposition 12. 1[ we can find a nested se- 
quence {Xj}^L Q G Fp, p > 2, each Xj having mesh norm hj := hxj satisfying 
jhj < hj + i < ^hj < ^j/i-o- If 4> G LP is an SBF, then define the Vj's to be 



These spaces have finite dimension equal to the cardinality of Xj and by virtue of the 
Xj's being nested, are themselves nested. At issue then are the Favard and Bernstein 
inequalities. Since any <fi that satisfies both will provide us with a Besov space via 
Proposition 16. 1 Jl we have the following result. 

Corollary 6.13 Let 1 < p < 00, <j)p = Gp + Gp * ip, where tp £ L 1 and < /3. Fix 
< 7 < (3 — n/p' and suppose that Vj = Q^„,Xj> with Xj as in \6.2(J\) . For allO < r < 
7 and all < r < 00, we have that f G B^. p if and only if {distx,p(/, Vj)}JL G b r r . 
The same conclusion holds true, with any 7 > 0, for all <f> that simultaneously satisfy 
\6. 3\) and i6.12]) . including the Gaussians, multiquadrics, etc. 

Proof: When Vj = Gcf>g,Xii t ne result follows immediately from the Bernstein in- 
equality in Theorem 16 . 1 1 and the Favard inequality in Theorem 16. 81 If 4> satisfies both 
(16.31) and (|6.12[) . then it also satisfies both the Bernstein inequality in Theorem 16.41 
and the Favard inequality in Theorem l6.10l As before, with the same set of V}'s, the 
same conclusion holds. □ 



Vj : c; c , v . j = 1,2 



, . . . , and V = {0}. 



(6.20) 
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6.4 Inverse theorems 



Inverse theorems give indications of rates of approximation being best, or nearly 
best, possible. We now establish inverse theorems for the approximation rates in the 
previous section and in [T5]. These involve Bessel-potential Sobolev spaces, and in 
addition Besov spaces. 

Theorem 6.14 Let 1 < p < oo and let (j) as in Theorem \6.1\ or Proposition \6.4\ If 

for f G L p , 1 < p < oo, or f G C(§"), p = oo, there are constants < < 7, ( € R, 
and Cf > such that 

dist LP(s „ ) (/,g , x )<c / (6.21) 
log 2 (^ x ) 

holds for all X G T p , then, for every < v < p,, f G If fOT]J holds for 

v = (i and some t > 1, then f G HP(S n ). Moreover, if in addition (f> satisfies the 
conditions in Corollary \6. 131 then for any r > > and < r < p, the function 
f is in the Besov space B T T p . 

Proof: Let the Vj's be as in (|6.20[) . and set fj := argmin (dist^s^ (/, Vj)), which 
always exits because Vj is finite dimensional. Since the Vj's are nested, we have 
that fj G Vk for all k > j. We want to show that fj is a Cauchy sequence in Hf,. 
From the Bernstein estimate in Theorem 16. II - or Proposition l6.4| - and the inequality 

hj+i/Qj+i < Pi wc have 

- fj\\ m < C P v hJ^\\fj^ - fih < Cp v hJ^{\\f j+1 - /||p + 11/ - fj\\ P ). 
And by (|6. 2 1[) . we also have 

\\f j+ i - fj\\m < CcfP V hJ+i( h ?+i iog 2 "*(^+i) + ^ logs *(*,-)) 

< Cc f p u h 2-^- v K j+ V ((h + j + l)- ( + 2»(h Q + j)-*) 

< C'c f 2-^-' /)j j- t 

where C is independent of j. Take k > j. Using the previous inequality and a 
standard telescoping-scries argument, we arrive at this: 

k 

||/,--/ fe || ff g<C"(^2-^) m m-*). 

m=j 

Letting j, k — ► 00, we see \\fj — fkWu^ ~ * when p > v and r G R or when p = v 
and t > 1 . Thus, fj is a Cauchy sequence in Hf, and is therefore convergent to 
/ G Hf,. Moreover, by (|6~2T|t with X = Xj, we sec that fj — > / in LP, so / = / 
almost everywhere. Hence, we have / G H%. The statement concerning Besov spaces 
follows from two things: the observation that aj := distil (/, Vj) < Cf2~Mj~ t , so 
||aj| rr < 00 whenever < r < p and rt > 1, and Corollary 16. 131 □ 

For the case v = p, < t < 1, the inverse theorem fails for Bessel-potential 
Sovolev spaces, but still remains valid for Besov spaces with r > t . 
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7 Concluding Remarks 



There are connections between this paper and [15[ 113] . In these papers quasi inter- 
polatory SBF networks were obtained yielding near best approximants for functions 
in Sobolev classes. The associated quasi-interpolation operators were constructed in 
the Fourier domain. The paper [15] focused on sequences corresponding to the c°° 
case treated within this paper. The paper [lOj dealt with sequences connected to the 
"perturbations of Green's functions" case. For example, let tp be a perturbation of 
a Green's function as described in this paper. If the Fourier coefheents of tp satisfy 
the "difference condition" as stated in , [T3j then it is in L 1 . The examples given in 
Section 3 satisfy both kinds of conditions. 

In [TS1 H3] , the quasi- intcrpolatory SBF networks were shown to give best results 
in the sense of n-widths. In this paper, using the frame approach, we have shown 
the quasi-interpolatory networks are also optimal for approximation of individual 
functions. Also note that in [13] . Marcinkiewicz-Zygmund measures generalizing the 
measure that associates fi q (R^) with each £ were introduced. These measures were 
used to derive Q21 Prop. 4.1 & (4.15)], which have overlap with the current Prop. 
4.. 4, Lemma 4.7 and estimate (5.4). 

In [lOj , the quasi-interpolation polynomial operators were further utilized to show 
that, in the presence of certain singularities, they exhibited better approximation 
properties than traditional methods. Also [TDl Prop. 4.3] is related to Proposition 
4.10 given here. Finally there is material closely connected to Theorem 4.1 appearing 
in [6, [HI Proposition 4.1]. Another version of the operator Bj was introduced in 

PI1 : a Af) = E?=o h i l / 2J ) P l(f)> where h : [°! °°) [°j 00 ) is a function in C k , equal 
to 1 on [0, 1/2] and on [l,oo). An early form of Theorem 4.1 was Theorem 3.4 of 
[M] . Frames, based on the crj(f) operator can be constructed as in jTTl |T8] using 
h(t) — h(2t) in place of k used in the construction given here. 

Finally we mention that the idea of using minimal separation for converse theo- 
rems and Bernstein inequalities goes back to [53], see also [TT]. Also, for the neural 
network community, we note that the number of neurons is not used as a measure of 
complexity, but rather the minimal separation of the nodes. 
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